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基于 遥感 与 气象 数据 的 冬小麦 主 产 区 籽粒 蛋白质 含 量 预报 
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摘 要 : 开展 小 麦 籽粒 和 蛋白质 含 量 的 监测 预报 研究 对 于 指导 农户 调 优 栽 培 、 企 业 分 类 收 储 、 期 货 小 麦 价格 、 


进口 政策 调整 等 具有 重要 意义 。 本 研究 以 冬小麦 主 产 区 (河南 省 、 山 东 省 、 河 北 省 、 安 徽 省 和 江苏 省 ) 为 
研究 区 域 ， 构 建 了 冬小麦 籽粒 蛋白 质 含量 多 层 线性 预测 模型 ， 并 实现 了 2019 年 冬小麦 蛋白 质 含量 预报 。 为 


了 解决 预测 模型 在 年 际 扩 展 和 空间 扩展 存在 偏差 的 问题 ， 在 蛋白 质 含量 估算 模型 中 考虑 了 气象 因素 (温度 、 
降水 、 辐 射 量 ) 、 冬 小 麦 筋 型 、 抽 穗 一 开花 期 增强 型 植被 指数 (EVI) 等 因素 。 结 果 表 明 ， 融 合 3 个 气象 因素 
的 蛋白 质 含量 估算 模型 建 模 集 精 度 (R = 0.39, RMSE = 1.04%) 与 验证 集 精 度 (R? = 0.43, RMSE = 0.94%) 
均 高 于 融合 2 个 气象 因子 的 估算 模型 和 单个 气象 因子 的 估算 模型 。 将 蛋白 质 含量 估算 模型 应 用 冬小麦 主 产 区 
的 蛋白 质 含量 遥感 估算 ， 得 到 了 2019 年 冬小麦 主 产 区 品质 预报 图 ， 并 形成 黄 淮海 地 区 冬小麦 品质 分 布 专题 
图 。 本 研究 结果 可 同时 为 后 续 小 麦 种 植 区 划 和 实现 绿色 、 高 产 、 优 质 、 高 效 粮食 生产 提供 数据 支撑 。 
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1 引言 


随 着 社会 经 济 的 日 益 发 展 和 生活 水 平 的 不 断 
提高 ， 当 前 人 们 对 农产品 的 需求 已 经 从 “有 没 
有 ”“ 饮 不 够 ”转向 “好 不 好 ”“ 优 不 优 ”。 但 是 ， 
中 国 农业 大 而 不 强 、 多 而 不 优 的 问题 依然 存在 。 
2019 年 ， 农 业 农村 部 、 国 家 发 展 和 改革 委员 会 、 
科学 技术 部 等 联合 发 布 的 《国家 质量 兴 农 战略 规 


收 稿 日 期 2021-03-22 ”修订 上 日期: 2021-04-25 


划 (2018—2022) 》 中 提出 “推进 农业 由 增产 导 
向 转向 提 质 导向 ”和 “ 走 质量 兴 农 之 路 ”等 发 展 
目标 。 籽 粒 蛋 白质 含量 是 小 麦 品质 评价 的 关键 指 
标 ， 而 品质 的 高 低 决 定 小 麦收 购 价格 、 加 工 用 途 
和 使 用 价值 等 O, Ke, FRENE MEAE 
量 的 监测 预报 研究 ， 对 于 指导 农户 调 优 栽 培 、 企 
业 分 类 收 储 、 期 货 小 麦 价 格 、 进 口 政策 调整 等 具 
有 重要 意义 “1。 
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基于 作物 籽粒 蛋白 质 含 量 
Content，GPC) 遥感 监测 预报 的 研究 目前 已 得 到 
探索 和 初步 应 用 ， 可 以 归纳 为 4 类 : (1) 基于 
“遥感 信息 -籽粒 蛋白 质 含 量 ” 模 式 的 经 验 模型 ， 
该 模型 通过 oe 
(敏感 波段 、 植 被 指数 、 红 边 参 数 等 光谱 特征 
直接 构建 作物 GPC 统 计 经 ee a EF 
“遥感 信息 -农学 参数 -籽粒 和 蛋白质 含量 ”模式 的 
定量 模型 ， oe Pe 
学 参数 之 间 的 定量 关系 及 农学 参数 与 GPC 之 间 
的 定量 关系 ， 构 建 GPC PM; (3) 基于 
遥感 数据 和 生态 因子 的 GPC 半 机 理 模 型 ， 该 模 
型 考虑 了 作物 氮 素 运转 机 理 和 生态 因子 对 GPC 
的 影响 ， 改 善 了 模型 的 年 际 扩展 性 和 空间 转移 
PESOS (4) 基于 遥感 信息 和 作物 模型 结合 的 机 
理解 释 模 型 ， 该 模型 综合 考虑 籽粒 和 蛋白质 形成 过 
程 中 各 种 生态 因子 的 影响 ， 通 过 遥感 信息 和 作物 
生长 模型 耦合 的 同化 方法 ， 调 整 模型 模拟 变量 与 
遥感 观测 值 的 误差 达到 最 小 ， 以 调整 作物 模型 的 
初始 参数 和 状态 变量 ， 进 而 实现 GPC 的 预 
测 "*。 综 合 分 析 以 上 各 类 方法 与 模型 的 选择 ， 
前 两 类 模型 研究 较 多 ， 操 作 简 单 且 易 实 现 ， 但 机 
理解 释 性 不 强 ， 在 区 域 间 和 年 际 间 扩展 应 用 过 程 
中 监测 结果 偏差 较 大 ; 同化 模型 方法 考虑 过 多 的 
输入 变量 ， 并 且 复 杂 的 同化 算法 、 品 质 生 长 模型 
本 地 化 精度 和 高 耗 时 运算 时 间 问 题 ， 限 制 了 该 类 
方法 的 区 域 大 面积 应 用 …。Li 等 和 Xu 等 六 
通过 综合 考虑 影响 品质 遥感 预报 模型 中 环境 变异 
因素 及 品种 筋 型 因子 ， 引 入 分 层 线性 模型 (Hier- 
archical Linear Models, HLM) 解释 籽粒 和 蛋白质 
含量 -遥感 -环境 的 般 套 问题 ,构建 基 于 开花 期 的 
冬小麦 品质 预测 模型 ， 很 好 地 解决 了 模型 年 际 扩 
展 和 空间 扩展 存在 较 大 偏差 的 问题 。 因 此 ， 发 展 
以 遥感 信息 〈 反 映 作物 养分 状况 ) 与 环境 因 
ieee 结合 的 GPC 半 机 理 预 测 方法 ， 
构建 普 适 性 及 适用 于 时 空 扩 展 的 籽粒 蛋白 质 遥 感 
预测 模型 具有 可 行 性 。 

本 研究 拟 在 前 期 构建 的 冬小麦 品质 分 层 遥 感 


(Grain Protein 


预测 的 基础 上 2” ， 初 步 尝 试 在 全 国 冬小麦 主 产 
区 (河南 省 、 山 东 省 、 河 北 省 、 安 徽 省 和 江苏 
省 ) 开展 GPC 模 型 的 构建 ， 以 期 实现 冬小麦 品 
质 收获 前 预测 。 


2 数据 与 方法 


2.1 冬小麦 研究 区 


试验 分 别 于 2008 年 、2009 年 和 2019 年 在 全 
国 冬小麦 5 个 主 产 省 (河南 省 、 山 东 省 、 河 北 
省 、 安 徽 省 和 江苏 省 ) 进行 数据 采集 ， 包 括 样 点 
取样 及 测试 、 小 麦 品 种 及 筋 型 确定 、 气 象 数据 收 
集 和 遥感 数据 获取 ， 共 获取 样品 898 个 ， 去 除 异 
第 点 后 得 到 样品 864 个 ， 其 中 2008 年 200 个， 
2009 年 283 个 , 2019 年 381 个 , 测 样 点 分 布 如 图 1 
所 示 。 其 中 ， 强 筋 小 麦 86 个 ， 强 筋 和 中 筋 混 合 
区 小 麦 249 个 ， 中 筋 小 麦 380 个 ， 中 筋 和 弱 筋 混 
合 区 小 麦 2 个 ， 弱 筋 小 麦 147 个 。 
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图 1 冬小麦 实测 样 点 分 布 图 


Fig. 1 Distribution map of the measured winter wheat 
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2.2 数据 获取 


2.21 冬小麦 籽粒 蛋白 质 测定 

在 冬小麦 成 熟 期 ， 获 取 农 户 冬小麦 样品 和 
样 点 经 纬度 ， 每 个 农户 收割 3 个 样 点 ， 每 个 样 点 
1m， 经 过 自然 蚊 晒 、 风 干 、 称 重 和 求 均 值 。 最 
后 利用 近 红 外 谷物 分 析 仪 FOSS Infratec™ 1241 
(Tecator， 赫 加 奈 斯 ， 瑞 典 ) 测定 冬小麦 GPC 
(14% FÆ) 并 记录 。 
2.2.2 ”冬小麦 品种 及 筋 型 确定 

在 冬小麦 样品 测试 过 程 中 ， 同 时 记录 农户 种 
植 冬小麦 的 品种 信息 ， 通 过 查询 品种 的 信息 确定 
品种 的 强 /中 / 弱 盘 特性。 区域 品质 预报 过 程 中 ， 
考虑 品种 的 多 样 性 及 复杂 性 ， 参 考 2001 年 农业 
部 发 布 的 《中 国 小 麦 品质 区 划 方 案 》 所 指定 的 依 
据 和 原则 ， 利 用 全 国 农 业 技术 推广 服务 中 心 收集 
汇总 的 品种 推广 数据 ， 以 及 以 县 为 单位 发 布 的 当 
年 主推 品种 ， 确 定 每 个 县 种 植 小 麦 筋 型 特性 。 最 
终 形成 以 县 为 单位 的 冬小麦 品质 区 划 图 。 
2.2.3 ”EOS/MODIS 还 感 影像 

以 抽穗 一 开花 期 为 时 间 节 点 进行 冬小麦 品质 
预报 ， 本 研究 选择 基于 中 分 状 率 成 像 光 谱 仪 
(Moderate-resolution Imaging Spectroradiometer, 
MODIS) 遥感 影像 构建 的 增强 型 植被 指数 (En- 
hanced Vegetation Index, EVI) 进行 建 模 呈 。 如 
公式 (1)。 

EVI =2.5x (Ryi—Rrea)/(Ryigt6*Rrei-7-5* Rene) (1) 

其 中 ，R\,、Rrw 和 Rai 分别 代表 近 红 外 波段 、 
红 光 波段 和 蓝光 波段 的 反射 率 。 为 消除 区 域 物 候 
差异 与 有 云 覆 盖 影 响 ， 分 别 获取 2019 年 和 2020 
年 4 月 一 5 月 MODIS 数据 ， 分 析 获 取 对 应 像 元 最 
大 值 作为 抽穗 一 开花 期 的 影像 结 
2.24 气象 数据 

冬小麦 研究 区 的 气象 栅 格 数据 从 欧洲 中 期 预 
报 中 心 (European Center for Medium-Range 
Weather Forecasts, ECMWF) (http: //www. ec- 
mwf.int/) 获取 。 本 研究 气象 数据 主要 包括 4 月 
份 的 月 积温 (°C + 日)、 月 太阳 辐射 量 (MJ/m?) 
和 月 总 降雨 量 (mm)， 选 择 这 3 个 气象 因子 进行 


冬小麦 品质 预报 是 因为 研究 区 内 冬小麦 从 起 身 期 
到 开花 期 约 为 一 个 月 时 间 〈 不 分 年 份 )， 而 该 生 
育 阶 段 包括 的 拔节 期 、 挑 旗 期 、 抽 穗 期 及 开花 期 
是 冬小麦 生长 最 重要 的 生育 时 期 ， 对 于 群体 的 建 
成 、 穗 粒 数 以 及 养分 积累 具有 决定 作用 。 利 用 
Python 3.7 (Python Software Foundation， 波 特 
=, 美国 ) 及 Matlab 2014 (Mathworks Inc， 纳 
Woe, FEL) 分 别 进行 数据 下 载 及 读 取 。 


2.3 冬小麦 品质 预报 模型 


分 层 线性 模型 (Hierarchical Linear Model, 
HLM) 是 一 种 考虑 了 数据 山 套 结构 (例如 学 生 
PERF BER) 的 最 小 二 乘 回归 分 析 ， 然 而 相 比 一 
般 最 小 二 乘 回归 分 析 ，HLM 在 针对 数据 之 间 的 
不 独立 性 时 ， 将 数据 集 分 层 、 综 合 分 析 层 内 数据 
(第 一 层 模型 ) 以 及 层 内 与 层 外 数据 (第 二 层 模 
型 ) 之 间 的 关系 "1。 

在 考虑 籽粒 品质 的 区 域 及 年 际 模型 构建 时 ， 
地 域 及 年 际 之 间 的 环境 因素 (第 二 层 模 型 ) 对 品 
质 的 影响 会 造成 遥感 信息 与 品质 模型 (第 一 层 模 
型 ) 的 差异 ， 因 此 考虑 HLM 模型 进行 品质 预测 
建 模 具有 可 行 性 “2 。 本 研究 中 ， 第 一 层 模型 是 
基于 GPC 与 植被 指数 和 将 粒 筋 型 构建 ， 公 式 
如 下 。 

GPC, = yt By EVI, +B, Glu +r; (2) 

其 中 ，CGPC ,为 籽粒 和 蛋白质 含量 ; Glu 为 冬 小 
麦 品种 筋 型 值 。 本 研究 依据 Li 等 ”的 研究 ， 强 
筋 、 中 强 筋 和 弱 筋 小 麦 品 质 分 别 取 值 1、2 和 3， 
后 期 在 区 域 应 用 过 程 中 有 强 筋 和 中 筋 种 植 区 与 中 
筋 和 弱 筋 种 植 区 ,分别 取 值 1.5 和 2.5。EVI 为 增 
强 型 植被 指数 ，Xu 等 ' 引 研究 表明 EV 下 与 冬小麦 
品质 较 其 他 植被 指数 具有 更 好 的 相关 性 ， 因 此 本 
研究 直接 选择 了 EV 指数 作为 模型 输入 参数 。 
Bo、PBi 和 pp; 分别 为 第 一 层 模型 的 截 距 、 植 被 指数 
系数 和 品种 筋 型 系数 。x 表 示 随 机 误差 。 第 二 层 
模型 基于 第 一 层 模 型 中 的 模型 系数 与 外 部 气象 数 
据 (温度 、 降 水 、 辐 射 量 ) 的 归 一 化 值 构建 ， 如 
公式 《(3)s 
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b,j= Yo + ynRad + ynTem + yanPre + u, (3) 
HP, nRad, nTem 和 npPre 分 别 为 每 个 栅 格 
4 月 份 月 太阳 辐射 量 归 一 化 值 、 月 积温 归 一 化 值 
和 月 总 降雨 量 归 一 化 值 ; 思 .分 别 对 应 第 一 层 模型 
中 截 距 、 增 强 型 植被 指数 系数 和 品种 筋 型 系数 ; 
yw 为 第 二 层级 的 截 距 ，y,,、ys 和 ys 分 别 为 第 二 层 
级 总 辐射 量 、 总 降水 量 和 积温 的 模型 系数 ，/ 心 则 
表示 为 随机 误差 。 气 象 数据 的 归 一 化 见 公式 (4) 
(以 Rad 为 例 )。 
nRad = (Rad- Rad „)/(Rad „— Rad œ) (4) 
Fp, Rad, Rad 和 Rad p 为 分 别 每 个 栅 格 
4 月 份 月 太阳 辐射 量 、 最 大 太阳 辐射 量 和 最 小 太 
阳 辐 射 量 ，MJ/m”。 为 进一步 测试 模型 ， 考虑 不 
同 的 气象 因子 组 合并 对 比 模型 的 精度 ， 以 确定 最 
优 的 气象 因子 。 模 型 的 构建 与 验证 通过 Matlab 
2014 实 现 。 


2.4 统计 分 析 


经 过 异常 点 剔除 后 获得 2008 年 、2009 年 和 
2019 年 样 点 数据 共计 864 组 ， 将 三 年 数据 混合 后 
随机 选择 其 中 80% 数据 用 于 建 模 〈 共 计 691 组 )， 
剩 下 20% 数 据 用 于 模型 精度 验证 (共计 173 组 )。 
利用 Matlab 2014 实 现 对 冬小麦 GPC 预测 值 与 实 
测 值 之 间 的 统计 分 析 。 采 用 决定 系数 RR 和 均 方 根 
误差 RMSE 验证 预测 模型 的 精确 度 与 可 靠 性 。R? 
值 越 接 近 于 1， 说 明 模 型 拟 合 程度 越 高 。RMSE 
值 越 小 ， 说 明 模 型 预测 能 力 越 强 ， 稳 定性 好 且 可 
靠 性 高 。 另 外 ， 为 进一步 分 析 冬 小 麦 品质 预测 模 
型 中 各 个 自 变 量 因子 对 最 终 GPC 的 贡献 度 ， 采 
用 相对 重要 性 分 析 法 (Dominance Analysis) "9 
对 模型 中 各 个 自 变 量 因子 对 GPC 的 相对 重要 性 
进行 计算 和 排序 。 


3 结果 与 分 析 


3.1 冬小麦 品质 区 划分 析 


参考 《中 国 小 麦 品质 区 划 方 案 》 所 指定 的 依 
据 和 原则 ， 结 合 5 个 冬小麦 主 产 省 各 目的 品质 区 
划 ， 以 及 当前 主推 品种 的 面积 及 市 县 分 布 ， 最 终 


形成 冬小麦 品质 区 划分 布 图 (图 2)。 
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图 2 2019 年 冬小麦 品质 区 划 专 题 图 
Fig. 2 Thematic map of winter wheat quality regionalization 


in 2019 


总 体 来 看 ， 强 筋 小 麦 主要 以 北部 地 区 种 植 为 
主 ， 如 河北 中 南部 、 河 南 省 北部 、 山 东 省 临沂 市 
和 滨州 市 等 。 中 筋 小 麦 以 山东 省 种 植 面积 最 广 、 
河南 南部 地 区 、 江 苏 省 北部 和 安徽 省 北部 等 也 有 
分 布 。 弱 筋 小 麦 则 主要 分 布 在 江苏 南部 及 沿海 地 
区 、 安 徽 中 南部 及 河南 南部 地 区 等 。 通 过 调查 发 
BR, 冬小麦 的 品质 区 划 除 了 与 生态 因子 (气象 及 
土壤 等 ) 相关 以 外 ， 也 与 品种 推广 和 政府 决策 相 
关 ， 例 如 2019 年 山东 省 济 麦 22、 和 鲁 原 502 ALLA 
28 号 等 中 筋 品 种 推广 占 总 面积 的 95.3%， FR PALL 
东 省 适合 强 筋 和 中 强 筋 小 麦 种 植 ， 但 综合 来 看 强 
筋 小 麦 品种 占 比 不 大 。 另 外 ， 在 政府 决策 方面 ， 
2019 年 河北 省 发 布 的 《河北 省 强 筋 小 麦 产业 提 质 
增 效 推进 方案 (2019 一 2022 年 ) 》 确 定 以 葛城 、 
WXEL. BRI. TORR. WAS. Bese. FEB. TE, 
南 和 、 沙 河 、 大 名 、 肥 乡 等 为 主要 推广 区 县 。 本 
研究 所 确定 的 小 麦 品质 区 划 专 题 图 是 在 小 麦 品质 
区 划 和 2019 年 各 市 县 的 主推 品种 调查 分 析 后 形 
成 的 。 在 确定 过 程 中 ， 受 县 级 单位 信息 发 布 不 
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全 、 部 分 市 县 品种 复杂 等 影响 ， 在 最 终 确定 品 
筋 型 过 程 中 ， 以 市 为 单位 确定 较 多 ， 后 续 的 研究 
及 调查 过 程 中 需要 更 精确 的 品种 推广 范围 ， 以 确 
定 准 确 的 冬小麦 筋 型 分 布 信息 ， 为 后 续 的 品质 预 
测 模 型 提供 可 靠 的 筋 型 空间 分 布 信息 。 


3.2 GPC 预报 模型 的 气象 因子 选择 


以 不 同 的 气象 因子 参与 分 层 线性 模型 进行 冬 
小 麦 GPC 预报 模型 构建 ， 并 对 比 模型 精度 进行 
气象 因子 的 评价 ( 表 1)。 结 果 表明 ， 单 独 一 个 气 
象 因子 进行 品质 预报 时 ， 以 ARad 为 输入 变量 构 
建 的 冬小麦 GPC 预报 模型 精度 (R = 0.31, 
RMSE = 1.12%) 均 高 于 分 别 以 n7Tem 和 7mPre 为 输 
入 变量 构建 的 冬小麦 GPC 预报 模型 。 两 个 气象 
因子 参与 冬小麦 GPC 预报 模型 的 精度 有 一 定 的 
提高 ， 以 nRad 和 nPre 组 合 构 建 的 冬小麦 GPC 预 
报 模型 较 好 (R = 0.38, RMSE = 1.05%)。 三 个 
气象 因子 共同 参与 冬小麦 GPC 预报 精度 最 佳 ， 
模型 的 RR 和 RMSE 分 别 为 0.39 和 1.04%。 考 虑 到 
冬小麦 籽粒 和 蛋白质 形 成 过 程 中 受 辐 映 量 、 温 度 和 
降水 的 综合 影响 ， 在 区 域 大 面积 模型 构建 过 程 
中 ， 气 象 因 子 存在 较 大 的 空间 变异 性 。 有 关 气 象 
因子 对 GPC 预报 模型 在 不 同 尺度 范围 上 的 影响 ， 
有 必要 在 品质 关联 的 气象 要 素 的 指标 筛选 及 关联 
度 分 析 等 方面 进行 进一步 研究 ， 特 别 是 气象 因素 
在 区 域 的 空间 变异 。 

进一步 通过 相对 重要 性 分 析 GPC 预报 模型 
表 1 不 同 气象 因子 构建 的 冬小麦 GPC 预报 模型 精度 
Table 1 Predicting GPC model in winter wheat by different 


meteorological data 


GPC 模 型 建 模 集 ”GPC 模 型 验证 集 


中 各 个 自 变 量 因子 对 GPC 预测 的 贡献 。 从 表 2 结 
果 中 可 以 看 出 ， 冬 小 麦 品种 对 GPC 的 影响 是 最 
大 的 ，Glu 相 对 重要 性 达到 75.31%， 品 种 的 遗传 
特性 直接 决定 品质 的 高 低 ， 考 虑 到 较 大 的 研究 区 
域 ， 品 种 的 多 样 性 对 模型 结果 具有 较 高 的 相关 
性 ; 气象 因子 对 GPC 的 重要 性 依次 为 辐射 量 
(nRad)、 温 度 (nTem) 和 降水 量 (nPre)， 相 对 
重要 性 分 别 为 10.05%、7.61% 和 6.54%; 遥感 信 
息 重 要 性 最 低 。 但 是 ， 当 以 局 部 区 域 (以 河北 省 
冬小麦 研究 区 为 例 ) 为 研究 区 时 ， 局 部 区 域 的 品 
种 多 样 性 减少 ，Glu 的 相对 重要 性 降低 到 
48.22%; 同时 由 于 局 部 区 域 的 气象 条 件 变 异 减 
小 ， 除 nRad 比重 相对 增加 外 ，n7Tem 和 nPre 依 次 
减少 ;反之 ,体现 作物 长 势 及 养分 差异 的 EVI 的 
相对 重要 性 增加 (22.65%)。 

表 2 冬小麦 GPC 预 测 模型 中 各 个 自 变 量 因 子 的 贡献 度 


Table 2 Contribution degree of each independent variable in 


GPC model of winter wheat 


模型 自 变量 黄 淮海 地 区 河北 省 
因子 相对 重要 性 /% 排名 相对 重要 性 /% 排名 
Glu 75.31 1 48.22 1 
nRad 10.05 2 20.32 3 
nTem 7.61 3 5.08 4 
nPre 6.54 4 3.73 5 
EVI 0.49 5 22.65 2 


通过 以 上 分 析 初 步 得 到 以 下 结论 。 模 型 中 各 
个 自 变 量 因 子 对 GPC 的 贡献 度 在 不 同 的 研究 太 
度 有 所 差异 : 局 部 区 域 品种 与 气象 因子 差异 较 
小 ， 植 株 长 势 及 养分 的 变异 对 GPC 的 结果 影响 
更 大 ; 反之 ,大 范围 区 域 品种 与 气象 因子 的 较 大 
差异 ， 导 致 植株 长 势 及 养分 的 变异 对 GPC HY 


编号 气象 因子 
R RMSE% Rè RMSE% 


影响 降低 。 当 然 ， 本 研究 在 样 点 收集 过 程 中 通 


1 nRad 0.31 1.12 0.37 0.99 
2 nlem 0.26 1.15 0.29 1.05 
3 nPre 0.31 1.11 0.28 1.06 
4 nRad,nTem 0.35 1.08 0.41 0.96 
5 nTem, nPre 0.33 1.10 0.29 1.05 
6 nRad, nPre 0.38 1.05 0.39 0.97 
7 nTem, nRad, nPre 0.39 1.04 0.43 0.94 


YE: nRad, nTem 和 nPre 分 别 为 4 月 份 月 积温 归 一 化 值 、 月 太阳 
辐射 量 归 一 化 值 和 月 总 降雨 量 归 一 化 值 


过 示范 区 (2008 年 和 2009 年 ) 和 农户 协助 收集 
(2019 年 ) 获得 ， 样 点 空间 代表 性 不 强 ， 对 结论 
的 解释 还 有 待 进一步 验证 。 另 外 ， 本 研究 气象 因 
子 选 择 4 月 份 数据 作 为 模型 输入 ， 主 要 考虑 到 该 
生育 阶段 包括 的 拔节 期 、 挑 旗 期 、 抽 穗 期 及 开花 
期 是 冬小麦 生长 最 重要 的 生育 时 期 ， 对 于 群体 的 
建成 、 穗 粒 数 以 及 养分 积累 具有 决定 作用 。 对 于 
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区 域 应 用 而 言 ， 由 于 区 域 间 该 生育 时 期 时 长 差 
异 ， 后 期 在 进一步 品质 预测 的 研究 中 ， 需 结合 区 
域 物 候 提 取 ， 进 而 精确 获取 冬小麦 关键 生育 时 期 
的 气象 数据 。 


3.3 冬小麦 GPC 预报 模型 


以 3 个 气象 因子 为 输入 变量 构建 冬小麦 GPC 
预报 模型 作为 初步 的 区 域 预 报 模 型 。 图 3 是 GPC 
估算 模型 建 模 集 和 验证 集 的 实测 值 与 估算 值 关 系 
Klo ÆRE RH 0.39, RMSE 1.04%, ATVI 
看 出 样 点 数据 中 高 佑 的 样 点 数据 多 于 被 低估 的 样 
点 数据 ， 尤 其 是 部 分 2009 年 北京 样 点 和 江苏 样 
点 ， 但 大 部 分 样 点 是 集中 在 1:1 线 的 周围 ; 验证 
集 的 尼 为 0.43，RMSE 为 0.94%; 所 有 样本 数据 
( 建 模 集 + 验 证 集 ) 实测 值 与 估算 值 的 RR 为 0.40， 
RMSE 为 1.03%。 结 果 表 明 ， 结 合 遥 感 与 气象 数 
据 的 分 层 线性 预测 模型 在 区 域 GPC 预测 上 具有 
可 行 性 ， 并 且 模 型 的 预测 稳定 性 较 高 。 


3.4 冬小麦 GPC 区 域 预报 


将 GPC 模型 由 样 点 数据 推 及 到 面 域 范围 ， 
得 到 2019 年 冬小麦 GPC 估算 图 (图 4)。 越 红 的 
地 区 代表 GPC 越 高 ， 绿 色 代 表 GPC 偏 低 ， 基 本 
呈现 从 北 到 南 冬 小 麦 品质 逐渐 降低 的 趋势 。2019 
年 黄 淮海 地 区 冬小麦 GPC 最 高 的 地 区 是 河北 省 、 
山东 省 北部 以 及 河南 省 中 部 和 北部 ， 基 本 在 
13.6% 以 上 ， 其 中 河北 省 大 部 分 及 河南 省 北部 
GPC 达 到 14.6% 以 上 ; 山东 省 西部 、 河 南 省 东部 
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图 3 GPC 预报 模型 估算 值 与 实测 值 关系 图 
Fig. 3 Relationships between the measured and estimated 
GPC from forecasting model 
模型 分 为 两 层 : 第 一 层 模型 是 基于 GPC、EVI 和 
冬小麦 籽粒 筋 型 构建 ; 第 二 层 模型 是 基于 第 一 层 
模型 中 的 模型 系数 与 外 部 气象 数据 (温度 nTem、 


与 西南 部 地 区 、 安 徽 省 北部 与 江苏 省 北部 地 区 
GPC SEAS ZE 12.3% ~13.6% Zi], GPC Ax f AY He 
区 是 江苏 省 东南 部 沿海 地 区 与 河南 省 东南 部 地 
区 ，GPC 基本 低 于 12.3%。 


4 结论 


本 人 研究 在 864 个 冬小麦 实测 品质 样 点 (2008 
年 、2009 年 和 2019 年 ) 基础 上 ， 利 用 分 层 线性 
模型 构建 了 冬小麦 GPC 估算 模型 对 2019 年 冬 小 
KEK (河南 省 、 山 东 省 、 河 北 省 、 安 徽 省 和 
江苏 省 ) 范围 内 的 冬小麦 GPC 进行 估算 。 佑 算 


降水 aPre、 辐 射 量 nRad) 的 归 一 化 值 构建 。 
GPC 估算 模型 建 模 集 精 度 尼 为 0.39，RMSE 为 
1.04%; 验证 集 精度 有 RR 为 0.43，RMSE 为 0.94%。 
将 GPC 佑 算 模 型 应 用 于 冬小麦 主 产 区 的 GPC 遥 
感 估算 ， 最 终 得 到 2019 年 冬小麦 主 产 区 品质 预 
报 图 。 结 果 证 明 ， 构 建 的 籽粒 蛋白 质 遥 感 预 测 模 
型 具有 较 好 的 普 适 性 和 时 空 扩展 性 。 

本 研究 以 MODIS 影像 (空间 分 辩 率 为 
250 m) 为 例 初步 在 大 尺度 黄 淮 海地 区 冬小麦 种 
植 区 开展 应 用 ， 粗 分 辩 率 的 和 遥感 影像 数据 在 冬 小 
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图 4 2019 年 冬小麦 主 产 区 GPC 预报 图 


Fig.4 GPC forecast map of main winter wheat producing 


areas in 2019 


麦 面积 提取 方面 存在 混合 像 元 的 问题 在 本 研究 中 
未 进行 深入 探讨 。 因 此 ， 在 接 下 来 的 研究 中 有 必 
要 开展 高 分 辨 京 影像 数据 (Landsat-TM, Senti- 
nel-2 或 者 高 分 系列 卫星 数据 等 ) 在 冬小麦 品质 
预报 的 应 用 。 
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Abstract: With the rapid development of economy and people's living standards, people's demands for crops have changed from 
quantity to quality. The rise and rapid development of remote sensing technology provides an effective method for crop monitor- 
ing. Accurately predicting wheat quality before harvest is highly desirable to optimize management for farmers, grading harvest 
and categorized storage for the enterprise, future trading price, and policy planning. In this research, the main producing areas of 
winter wheat (Henan, Shandong, Hebei, Anhui and Jiangsu provinces) were chosed as the research areas, with collected 898 
samples of winter wheat over growing seasons of 2008, 2009 and 2019. A Hierarchical Linear model (HLM) for estimating 
grain protein content (GPC) of winter wheat at heading-flowering stage was constructed to estimate the GPC of winter wheat in 
2019 by using meteorological factors, remote sensing imagery and gluten type of winter wheat, where remote sensing data and 
gluten type were input variables at the first level of HLM and the meteorological data was used as the second level of HLM. To 
solve the problem of deviation in interannual and spatial expansion of GPC estimation model, maximum values of Enhanced 
Vegetation Index (EVI) from April to May calculated by moderate-resolution-imaging spectroradiometer were computed to rep- 
resent the crop growth status and used in the GPC estimation model. Critical meteorological factors (temperature, precipitation, 
radiation) and their combinations for GPS estimation were compared and the best estimation model was used in this study. The 
results showed that the accuracy of GPC considering three meteorological factors performed higher accuracy (Calibrated set: 
R? = 0.39, RMSE = 1.04%; Verification set: R? = 0.43, RMSE = 0.94%) than the others GPC model with two meteorological fac- 
tors or single meteorological factor. Therefore, three meteorological factors were used as input variables to build a winter wheat 
GPC forecast model for the regional winter wheat GPC forecast in this research. The GPC estimation model was applied to the 
GPC remote sensing estimation of the main winter wheat-producing areas, and the GPC prediction map of the main winter 
wheat producing areas in 2019 was obtained, which could obtain the distribution of winter wheat quality in the Huang-Huai-Hai 
region. The results of this study could provide data support for subsequent wheat planting regionalization to achieve green, high- 
yield, high-quality and efficient grain production. 


Key words: winter wheat; grain protein content (GPC); remote sensing; hierarchical linear model (HLM); meteorological data 
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